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ABSTRACT 


A  method  for  the  analysis  of  a  transonic  flow  field  in 
a  nearly  circular  duct  with  gradual  changes  in  crossi-s action 
is  developed.  A  linearized  differential  equation  for  the 
deviation  from  the  results  of  the  one-dimensional  theory  is 
presented  and  is  solved  by  the  superposition  of  particular 
solutions  obtained  by  a  product  hypothesis.  The  potential 
equation,  in  this  case,  was  simplified  for  the  vicinity  of 
the  sonic  velocity. 

Application  to  a  circular  subsonic-supersonic  duct  is 
made  in  order  to  determine  the  magnitude  of  propagated  sub¬ 
sonic  disturbances  and  the  resulting  velocity  distribution. 
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SECTION  I 


INTRODUCTION 

This  report  presents  an  analytical  method  for  the  calculation  of 
the  flow  field  in  a  nearly  circular  duct  with  relatively  gradual  changes 
in  cross-section.  Due  to  simplifications  introduced  by  means  of  the 
transonic  law  of  similarity ,  the  method  in  its  present  form  would  be 
restricted  to  flow  fields  with  Mach  numbers  close  to  unity. 

The  investigation  was  initiated  in  connection  with  a  plan  for 
modifying  the  10- Foot  Wind  Tunnel  at  Wright-Pat terson  Air  Force  Base. 

It  was  intended  to  vary  the  cross-sectional' area  with  small  inserts, 
rather  than  by  a  variation  of  the  radius.  A  method  was  evolved  to 
obtain  the  magnitude  of  deviations  from  the  ideal  parallel  flow  which 
might  be  caused  by  such  an  insert  design.  Before  the  actual  computa¬ 
tions  were  made  the  design  was  modified.  However,  the  method  itself 
is  of  interest,  in  that  there  are  many  other  applications.  An  appli¬ 
cation,  referring  mainly  to  a  modified  test  section  design,  will  be 
given  in  the  second  part  of  this  report. 

The  present  analysis  begins  with  a  basic  flow  as  given  by  one¬ 
dimensional  theory.  Any  deviations  from  the  basic  flow  are  considered 
small,  so  that  they  can  be  determined  by  a  linearized  partial  differen¬ 
tial  equation.  The  general  solution  to  the  differential  equation  is 
expressed  by  the  superposition  of  particular  solutions,  each  of  which 
is  the  product  of  a  Bessel  function  dependent  on  the  radial  distance 
y  from  the  duct  axis,  a  trigonometric  function  of  the  angular  position 
<o,  and  a  function  of  the  axial  position,  x.  Of  special  interest  is 
the  last  function  which  is  influenced  by  the  duct  shape. 

If  one  is  interested  only  in  supersonic  flow  analysis,  an  alternate 
method  of  computation  is  the  well-known  method  of  characteristics.  The 
method  Of  characteristics  is  more  general,  in  that  the  condition  of  a 
slow  change  in  cross-section  is  not  required.  The  method  of  character¬ 
istics  for  three  dimensions  is  similar  to  the  present  method  in  the 
requirement  of  a  nearly  circular  cross-section.  The  major  advantages 
of  the  present  method  are  that  the  subsonic  influences  can  be  deter¬ 
mined  and  that  being  essentially  more  analytical,  a  better  picture  of 
the  flow  behavior  in  general  cases  may  be  obtained.  In  some  cases, 
depending  upon  the  duct  shape,  the  solutions  can  be  obtained  in  a 
less  tedious  manner. 

Briefly,  the  report  starts  with  the  simplification  of  the  general 
flow  differential  equation  by  means  of  the  transonic  law  of  similarity. 
Then  a  general  method  of  solution  is  developed,  and  finally  this  result 
is  applied  to  some  specific  example. 
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SECTION  II 


THEORETICAL  CONSI DBR.1TI ONS 

Simplification  of  the  Plow  Differential  Equation  Near  the  Sonic  Speed 

This  analysis  is  based  on  a  cylindrical  coordinate  system  (Figure  l). 
Let  the  x  axis  be  the  duct  center  line,  y,  the  radial  distance  from  the 
x  axis,  and  u>  ,  the  angular  position  measured  clockwise  from  the  vertical 
position.  The  velocity  potential  is  denoted  by  $  •  Partial  derivatives 

are  denoted  by  subscripts  (i.e.,  M-r  (be) 

4  X  '  * 

The  general  three-dimensional  equation  for  a  steady  potential  flow 
(derived  in  Appendix  A)  is  given  by 

(1.D  <*«(«-£)+ $  ♦  Vao 

If  deviations  from  a  parallel  sonic  flow  are  small,  the  simplified 
transonic  equation  results.  (Derivation  may  be  found  in  Appendix  B).  • 

(1.2)  -U+0<P*x  <P*  +  <PYy  +  ^  ^  =  0 

where  ^  -  jL 

(1-5)  Vk.fc&rit-SXY 

and  r0  is  the  radius  at  the  throat. 

The  Pressure  and  Density  Relations 

The  new  potential  term  CP(x,y)  is  directly  related  to  the  velocity. 

However,  one  may  be  interested  in  the  pressure  and  density  relations  at 
various  positions  along  the  duct.  Bernoulli's  law  for  compressible  non- 
viscous  gases  is  given  by 

(1  .h)  ^  +  U4U  +  VAV  =  0 

The  last  term  is  negligible,  based  on  earlier  assumptions. 

Because  Au  =  <9*  c?  ,  one  can  write 

*  * 

Since  -  Q-  ,  the  preceding  equation  is  rewritten  as 
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(1.5) 


The  Boundary  Conditions 

Let  7  be  the  value  of  y  at  a  point  on  the  duct  surface  so  that 
r  *  r  (x,<*>).  This  function,  in  the  case  of  a  vertical  line  of  symmetry 
is  expressed  as 

(l  .6)  r  -  00  +•  £  (f*  (x)  cos  *  to) 


The  functions  r0,  r^,  T2, . rn  are  obtained  by  a  Fourier  analysis 

at  the  specific  x  positions,  and  are  written  as 


(1.6a)  ro (x)  =  Jz. 


5 


XXX 

?(%,<»>)  do> 

O 


(1.6b)  r* 


cos  vyo)  d  to 


Let  r  *  r  . 
ro 

The  boundary  condition  of  zero  flow  normal  to  the  duct  surface  must 
be  fulfilled.  If  the  differences  from  the  sonic  velocity  are  small, 
this  condition  is  expressed  as 

(1.7)  -  ^*<PY=0 

(The  derivation  may  be  found  in  Appendix  C). 

Let  us  introduce  for  the  slope 

(1.6)  -©•=£; 

Then  one  has  the  relation 

(1.9)  <Py='©o  +  E*ey  COS  VCO 

ii/here  ■©-v  is  an  additional  slope  for  the  unsymmetrical  configurations . 

The  Basic  Hypothesis  for  the  Solution 

The  solution  is  based  on  one- dimensional  flow  theory,  in  which  the 
velocity  distribution  over  each  cross  section  is  considered  constant. 

The  cross  section  is  written  as  functions  of  the  distance  from  the 
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minimum  cross-section.  The  potential  of  this  basic  flow  is  denoted  by1  <Pe 

,  where  ^*=q>e.  =  f<r,  yj').  Specifically,  as  derived  in  Appendix 
D,  this  gives 

»*“> 

The  entire  potential  is  expressed  by 

(1.11)  Cpr  CPe^+ 

The  CP  represents  the  deviations  from  the  basic  flow,  and  shall  be 
considered  small  enough  that  higher  order  terms  can  be  neglected.  If  the 
values  of  the  complete  potential  (Equation  (1.11))  is  inserted  into  the 
transonic  differential  equation  (Equation  (1.2))  and  second  order  terms 
of  are  neglected,  the  following  inhomogeneous  equation  is  given* 


For  the  representation  of  <P  ,  let  us  introduce  the  following  functions. 
Let  a  Bessel  function  of  V  ^  order  and  of  the  argument  be  denoted  by 
Let  be  the  argument  (where  ■  1,  2,  3  .  .  .)of  the  Bessel 

function  J%.  ,  for  which  ii  o  .  The  functions  2vvkare  given  by 

avj  ” 

(1.13)  Jv(y  7}v>0 

The  functions  2 Vy.  fulfill  the  following  differential  equation 

(i.ii,)  +C'5jV-vs)£v)‘=o  7  - 

There  exists  a  relationship  of  orthogonality  between  the  ?Vu.  >  80 
that  for  )k  v 

(1.15)  J  Y  <iy  =0 

Let  us  introduce 

(1.16)  [  Y  <iy-Kv>x 

^0 

A  hypothesis  for  the  perturbation  velocity  potential  is  chosen  in  such  a  way 
that  the  boundary  conditions  are  fulfilled. 
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(1.17) 


+  S(^YV-©'v(x^  cos  vco) 

+  (<V(x>  2orCY>)  +  to  ^v^(y)  cos  vw) 

The  functions  and  a.-yjj.C'iO  (  V  ■  1,  2,  3  •  •  •  )  must 

be  determined  in  such  a  way  that  the  differential  equations  for  (Q 
(Equation  (1.12)  is  fulfilled.  Inserting  this  hypothesis  into  Equation 
(1.12),  one  obtains 


(1.1S) 


-(*.»>[ <p„,  I( ( °-v  av)  +  <(>.„ i, (°-l 

+  a©-.  £vr  coi  vu')+  ^b.. 

I,  (*'r  Coi  v“>]  *J  V>-  (  ^  '»*  «•» 

♦  -  £  Zv^  c.s  »«)  -  (t* -)  [ T*  *  y” 

cos  vu»(cpBx  ■©‘v  +•  «C))=  -(*+»}(<*•.•<&«) 


All  the  terms  of  this  equation  which  contain  trigonometric  terms  of  j/'u 
(where  V  *  1,  2,  3,  U,  etc.)  must  vanish  separately. 

For  V  *  0  one  has  the  condition 

_  *•  *  •* 


A*>*» 


(1,19)  + 1  ^  v  c 

=.-z-e-.- Os.t>  S’.,  cpe„ 

The  conditions  for  V  ■  1,  2,  3,  .  .  .  are 

a.,V  <?*„  T, 


(1.20) 


*!,  «•*.(  ^ *  f  4**  -  T-  ^ 
=  05*) 

It  can  readily  be  seen  from  Equation  (1.1Q)  that 

(i.2i)  ee0=  cpB^  cpBxll 
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Further  simplifications  are  made  using  Equation  (1.11+)  with  Equa¬ 
tions  (1.19)  and  (1.20),  respectively. 

(1.22)  J, C <•<?».  +  ®B.t  *■-*)] 

-  £  1 X  ^  <5  •-  ^  (®..  <  *  «.«  «L 1 

(1.23)  -(*-h)  Jt  [fcvr(q>Bil  a;r  +  q>B^ av^] 

-  [V  a*/»]  ”  1  C^Pbx  ■Q'v) 

The  equations  for  the  perturbation  functions  are  obtained 

by  multiplying  with  the  term  y  and  integrating  between  y  *  0  to 

y  ■  1.  .Because  of  the  orthogonality  relations  for  the  Bessel  functions, 
all  witii  the  terms  J*  y  2.v^2.y/,  <Jy  drop  out  except  for  >*=V  ♦ 
Inserting  these  relations  in  Equations  (1.22)  and  (1.23)  one  obtains 
Equations  (1.2i|.)  and  (1.25),  respectively. 

°V  *  a-r  +  %-  °->0 

=  <P*«'©0 

K„^-  ("a.^  ys-iy  =  h  X(v  %.) 

'J0  |or 

(Derivation  may  be  found  in  the  Appendix  E).  Also, 

<1-27)  [  cpB,  a,;  +  cp8„  aY;  +  V  a,,.] 

(1.28)  X.,,  =  j"  e,r  (y)  y  <>y = "C Ov  V) 

(1.29)  £vr(y')  y'H”  (%.  fl 

From  Equations  (1.21+)  and  (1.27),  and  the  boundary  conditions  which 
will  be  given  later,  the  values  of  Ot^can  be  found.  The  resulting 
expression  for  the  perturbation  potential  (Equation  (1.17))  is  written  in 
a  more  convenient  manner  if  one  expresses  the  1/2  y^  and  ^  yv  by  means 
of  2;,,^  and  2V«.  .  ( 

"k  Y v  -  %  f:  ^  0 


(1.21+) 

where 

(1.25) 

(1.26) 
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where 

i  ^ 

°°r-  "  Z  y 

Let 

us  introduce 

(1.30) 

i  * 

-t*  -  °-oy~ 

(1.31) 

bvj*  t*  J* 

The  perturbation  potential  is  then  based  entirely  in  terms  of  the 
functions  2.  v  so  that 


(1.32) 


*  j|0  Ca*r  *JL  ?=.  t 


One  can  then  write  the  differential  equations  for  by  the 

introduction  of  Equations  (1.30)  and  (1.31)  into  Equations  (1 .24)  and 
(1.27). 

(1.33)  4>„ .  +  cpe„  ^  a! 


(1.34)  G'V)* 


-%r  // 

CX,lvi  *<P, 


♦  / 

<or 

♦  / 


=  “in  ike-. 


2Xor 


-  X 


f\  4-  a*  =  *yy*  ^vV 

® **  °“vr  +  'S+l  <X'T‘  'Ur\  ^ 


cp  and  its  first  derivatives  with  respect  to  x  are  continuous. 
Therefore,  according  to  this  relation,  snd  their  first  derivatives 

are  continuous. 

One  might  ask  why  the  hypothesis  (Equation  (1.32))  was  not  used 
from  the  beginning  instead  of  Equation  (1.17).  The  reason  is  that  the 
second  derivative  with  respect  to  y  does  not  converge  (because  of  the 
character  of  the  function  yv  ).  Therefore,  all  of  the  considerations 
regarding  the  fulfillment  of  the  differential  equation  for  Cp  would 
have  been  meaningless . 

Method  of  Solution 

# 

The  solution  of  the  perturbation  functions  for  the  differen¬ 

tial  equations  (1.33)  and  (1.34)  is  found  for  ary  duct,  within  the  defined 
limitations,  by  means  of  numerical  integration.  For  suitable  configurations, 
analytical  solutions  can  be  found.  In  some  cases,  a  piecewise  analytical 
solution  may  be  satisfactory.  In  any  case,  the  pressure  distribution 
along  the  axis  is  influenced  only  by  the  expressions  with  V  ■  0,  since 
the  nature  of  higher  order  Bessel  functions  is  such  that  there  is  no 
effect  along  the  axis. 

Deviation  from  a  parallel  flow  are  found  from  Cp^  .  The  only  expres¬ 
sions  which  are  of  Importance  along  the  axis  are  those  for  V  ■  0  and 
V  -  2.  Actually,  an  infinite  number  of  terms  are  needed  to  describe 
the  flow  completely,  but  the  higher  terms  in  y-  and  v  are  very  small 
and  are  of  influence  only  close  to  the  wall. 
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Conditions  for  the  Solution  of  Equations  (1.33)  and  (1.34) 


Equations  (1.53)  and  (1.34)  are  of  the  second  order,  and  consequently 
the  solutions  must  be  determined  by  two  boundary  conditions.  The  first 
condition  is  obvious.  With  the  assumption  of  a  subsonic  duct  extending 
to  negative  infinity,  one  may  assume  that  the  values  of  the  perturbation 
terms  should  vanish  as  one  approaches  negative  infinity. 

The  character  of  the  second  condition  is  less  obvious.  One  may 
first  ask  at  which  point  such  a  condition  might  be  given.  Certainly  it 
cannot  be  a  point  in  the  supersonic  region,  for  the  fulfillment  of  the 
condition  will  affect  the  entire  solution,  including  points  upstream  of 
the  point  where  the  boundary  condition  is  given.  Since  this  would 
violate  the  law  of  the  forbidden  signals,  the  location  of  the  other 
condition  should  be  expected  to  be  in  the  subsonic  region.  Actually, 
mathematical  reasoning  leads  to  the  location  of  the  boundary  condition 
at  the  narrowest  cross-section.  If  one  rewrites  Equations  (1.33)  and 
(1.34)  in  such  a  way  that  the  coefficient  of  the  highest  derivative  of 
the  term  atL.  is  one,  the  coefficients  of  the  other  terms  become  singular 
at  the  sonic  velocity.  The  solution  of  the  differential  equation  may  be 
singular  at  any  point  where  these  coefficients  are  singular. 

Consider  now  only  the  solutions  of  the  homogeneous  part  of  the  dif¬ 
ferential  equation,  since  the  influence  of  change  in  the  arbitrary 
constants  and  consequently,  the  influence  of  the  boundary  condition  is 
limited  to  that  part.  Let  the  basic  flow  velocity  be  represented  by 

** .  Then  the  two  linearly  independent  solutions  of  the  homogeneous 
differential  equation  are  represented  by 

(1.55)  -  P  (>) 

(1.36)  oXr  -  POO 

where  P  denotes  a  power  series  of  x,  which  starts  with  the  absolute  term. 

As  the  term  approaches  zero  (at  the  throat)  we  have  «  >  °  .  Then 
the  expression  ,  from  Equation  (1.36)  will  tend  to  infinity  as  '*•-*■0  . 

Since  this  is  physically  impossible,  the  second  condition  is  that 
remains  finite  at  the  throat. 

To  sunmarize  the  boundary  conditions,  the  flow  deviation  expressions 

Mb 

must  be  finite  at  x  *  0  and  zero  as  ®°. 

SECTION  in 
APPLICATIONS 

Flow  Through  A  Circular  Subsonic-Supersonic  IXict 

We  apply  this  theory  to  an  axially-symraetric  transonic  duct  configura¬ 
tion.  This  duct  is  similar  to  the  modified  10-Foot  Wind  Tunnel  test  section 
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at  Wright-Patt arson  Air  Force  Base.  The  duct  is  first  cylindrical 
(Section  A,  Figures  h  and  5)  followed  by  a  throat  region  of  parabolic 
profile  (Section  B),  and  finally  the  duct  expands  conically  (Section  C). 
The  configuration  is  geometrically  determined  by  the  subsonic  junction 
position  X.  ,  the  radius  of  curvature  R,  and  the  angle  of  expansion 
It  is  desired  to  trace  the  magnitudes  of  the  x-velocity  deviations 
from  a  uniform,  parallel  flow  through  the  duct  channel . 

In  the  differential  equation  (1.33)  for  ,  the  shape  of  the 

duct  expresses  itself  by  the  function  W  ,  which  occurs  in  the 
coefficients  of  a.* ''  and  ,  and  furthermore,  in  the  function 

which  occurs  in  the  righthand  term. 

In  Section  A,  the  duct  is  cylindrical  and  consequently  the  radius, 
slope,  and  terms  are 

(2.1)  r=rv-l+ 

(the  value  of  can  be  seen  from  the  subsequent  formula  (2.U)) . 

(2.2)  -©;=  O 

(2.3)  (for  section  A) 


A'he  profile  of  Section  B  is  an  approximation  for  a  circular  arc 
with  the  radius  of  curvature  R.  When  approximated  by  a  parabola  with 
the  vertex  at  the  minimum  section  of  the  duct,  the  shape  is  expressed 


by* 

(2.h) 

r  -  '  *  Tr- 

(2.5) 

'e’“  “  ^ 

(2.6) 

(for  section  B) 

Section  C  starts  tangential 
conically  at  the  angle  ©<  .  The 

to  the  parabolic  Section  B  and  expands 
shape  terms  are 

(2.7) 

(2.8) 

Oo  -  °< 

(2.9) 

© 

* 

n 

s 

it. 

(for  section  C) 

These  values  inserted  into  the  differential  equation  (1*2S)  lead 
to  the  specific  differential  equation  for  each  section.  v 


Section  A 

(2.10)  CpB< 


♦  n 
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Section  B 

Jfr  "  i 

(2.11)  *  0-or  -f  Ck0r  -v-  E,r  =  t 

The  constants  are:  _  x 

7)  »»* 

r*  \+*i 

(2.12)  ^  v»„  v  t 

W-*-l  R  ) 

(o  13)  C  —  Ha.  y-H  x«^_ 

^2.15;  -  aR  /_£_  r-vt 

\»+»  R/ 


Section  C 

»•"*>  [&,  <*  -  ??.)]*  <c + 1  (  £t)V 

The  solutions  of  the  differential  equations  (2.10),  (2.11),  and 
(2. lit)  will  now  be  determined.  The  complete  solution,  in  each  case, 
is  the  sum  of  the  general  solution  of  the  homogeneous  portion  and  a 
particular  solution  of  the  non-homogeneous  equation.  The  constants 
will  be  determined  later. 

Section  A 

The  solution  of  Equation  (2.10)  is  „  z  x 

/_  y  *  -/-  i* — v  * 

,  .  *  _r  cpex(vo;  .  r 

(2.15)  +  v-ij*  e 

Section  B 

The  particular  solution  of  Equation  (2.11)  is 

(a.16)  <x^  =  E,r+  E<rx 

where  x 

..  (-L-  is-y- 

(2-“> 

The  homogeneous  part  of  Equation  (2.11),  through  a  simple  trans¬ 
formation,  is  of  the  form 

(2.19)  y^y’ -[*+&)*]  y=° 
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the  solution  of  which  is 
(2.20)  y  =  Zp[2AX*] 

The  term  denotes  any  linear  combination  of  Bessel  functions  of 
order  p.  The  general  solution  of  the  homogeneous  part  of  Equation  (2.11) 
then  is  such  that  the  complete  solution  becomes 

(a.si)  =  EvE^*  +  C3r3.|>  Br-^]*CVN.[2^*4] 

Section  C 

A  particular  solution  of  Equation  (2.II4)  in  Section  C  is  given  by 


(2.22)  O,  -  P-. 


°A*  II., 


In  order  to  find  the  general  solution  of  the  homogeneous  part  of 
Equation  (1.28),  the  transformations  Us  (<-  and  Os^(*)=c{^u) 

are  carried  out.  The  differential  equation  is  then  written 

(2.25) 


where 


(2.2 1*)  £, 


1)3  t 


Equation  (2.23)  has  the  form 

(2.25)  V  bxy  -  o 

where  the  solution^  is  given  as 

(2.26)  V  =  £x  [t  V>" 

Consequently,  the  solution  of  Equation  (2.23)  is 

(2.27)  C|  -U'  £xt,3-E5y*U\] 


1  Jahnke ,  E.  and  Emde,  F.,  Tables  of  Functions .  Fourth  Edition.  Dover 
Publications,  19^5*  P»  11*7. 

2  Ibid.,  p.  li*7. 
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C9fK~  +  ax!z*l~fr  ~zQ  J^fe,)[g4M~Xfe»)  *>^J 

where  (*  »~  avo)4  ES/t*[ T^Ce.)  J^.(e,)  +J_±(e,)  L^.(e,)J 

P  -  9  c:  T  \/  4r 


eo=2E./V 

By  the  use  of  the  previous  formulae,  the  expressions  is 

determined  completely.  The  quantities  X«  ,  R  ,  and  which 

determine  the  duct  geometry  are  determined  by  the  specific  problem 
under  consideration. 


Discussion  of  Results 

In  the  previous  formulae  the  quantities  "*•  ,  ^  ,  and  **  are 

written  in  the  general  form.  This,  of  course,  enables  us  to  determine 
the  solutions  for  different  configurations. 

In  the  following  analysis  the  values  which  correspond  to  the  Wright- 
Batterson  Air  Force  Base  10— Pbot  Wind  Tunnel  will  be  inserted.  Before 
this  is  done,  the  effect  of  the  choice  of  the  subsonic  junction  position 
X,  is  studied.  Most  investigations  of  a  supersonic  flow  are  made 
disregarding  any  influences  of  the  subsonic  portion  of  the  duct,  since 
there  is  no  theoretical  basis  for  its  consideration.  It  is  of  interest 
to  determine  some  numerical  values  of  the  subsonic  effects.  Therefore, 
a  flow  pattern  determination  will  be  made  of  a  configuration  which 
corresponds  to  the  supersonic  part  of  the  10-Foot  Wind  Tunnel,  and 
three  different  subsonic  junction  positions. 


Accordingly,  the  duct  shape  gives 

R  -  50  r0  o<  -  .15  degrees 

Xq  »  -0.10,  -0.20,  -0.^0 


The  values  of  'A 

I 

DwtyS*0and  are 


are  the  Bessel  function  arguments  for  which 


-  3.S32 

V-  7.015 
^,-10.1735 
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The  values  of  >  and  ^  +f+  are  determined  ty 

Equations  (2.1 2),  (2.13),  (2.17),  and  (2.1S),  respectively.  The  junction 
constants  ,  Ciy^  ,  and  C6^  are  obtained  by  means  of  Equations 

(2.33),  (2.34),  and  (2.35),  respectively. 

Of  special  interest  is  the  velocity  distribution  along  the  axis  of 
the  duct.  The  x-velocity  terns  found  from  one-dimensional  theory  are 
given  try  Equations  (2.3),  (2.6)  and  (2.9).  The  additional  x-velocity 
along  the  axis  caused  by  the  perturbation  potential  are  found  from 
Equation  (1.32).  Since  3%(«) = I ,  one  obtains 

(2.36)  ^ 

#  ' 

Figure  6  shows  the  °~®,  values  through  Section  B  for  X„r 
In  spite  of  appreciable  initial  flow  deviations  at  some  distance  upstream 
from  the  throat,  the  effects  are  hardly  noticeable  at  the  minimum  cross 
section  and  in  the  supersonic  region.  A  corresponding  behavior  is  found 
for  the  terms  o.*^,  a*  *  ,  etc.  Accordingly,  even  extreme  changes  in 
the  subsonic  region  are  without  importance  for  the  supersonic  region. 

Figure  g  shows  the  velocity  distribution  along  the  axis  according 
to  the  one-dimensional  theory,  and  a  velocity  which  includes  the  correc¬ 
tions  given  by  ,  and  <*?, .  Beyond  x  -  0.5  the  velocity  distri¬ 

bution  is  rather  smooth,  showing  the  present  configuration  to  be  satis¬ 
factory.  The  ntfttrimnm  velocity  deviations  from  the  one-dimensional 
theory  have  a  magnitude  corresponding  to  a  Mach  number  of  about  .004. 

The  velocity  deviation  at  about  x  •  0.4  is  rather  large,  but  this  is  to 
be  expected  from  the  following  reasoning.  Point  M  (Figure  3)  is  con¬ 
sidered  to  be  the  junction  between  Sections  B  and  C.  The  effect  of  the 
change  of  shape  at  this  point  propagates  to  the  axi3  along  the  Mach 
lines  UN.  Upstream  of  MN  the  flow,  especially  along  the  axis,  corresponds 
to  the  parabolic  surfaces.  Consequently,  any  change  in  pressure  distri¬ 
bution  caused  by  the  transition  is  felt  only  downstream  of  point  B. 

The  one-dimensional  theory  (Figure  g,  curve  )  shows  the  influence  of 

the  junctions  to  be  right  at  point  A.  The  superposition  of  the 
value  shows  the  occurrence  of  the  change  at  x  0  .38.  The  x-coordinate 
of  Point  N  in  Figure  3,  oomputed  by  an  exact  relationship  for  the  Mach 
angle,  is  x  •  .367.  Figure  g,  curve  shows  ’’hat  might  be 

expected  from  a  more  general  reasoning.  The  functions  <^*, ',  , 

and  a*,'  are  given  in  Figure  7. 

A  Remark  on  the  Application  to  a  Non-Symmetric  Configuration 

Although  no  examples  of  the  application  for  an  unsymmetric  flow  are 
carried  out,  an  idea  of  the  procedure  will  be  given.  The  shape  of  the 
duct  must  be  expressed  according  to  liquation  (1.6).  For  instance,  one 
must  carry  out  a  harmonic  analysis  of  the  function  r(x,»)  with  respect 
to  «0  for  various  values  of  x.  This  can  be  done  by  a  numerical  or 
analytical  evaluation  of  the  integrals  (1.6a)  and  (1.6b),  or  by  means 
of  a  harmonic  analyzer  .  Then  the  values  of  are  found  from 

Equation  (1.9).  Following  this  preparation,  the  procedure  can  be  carried 
out  in  the  same  way  as  for  an  axially  symmetric  duct. 


Solution  with  as  A  Function  of  X.* 

A  solution  of  the  differential  equation  is  determined  for  the 
coefficient  as  a  power  of  the  longitudinal  position.  The  solution 

and  its  derivation  will  be  found  in  the  Appendix  (Section  IV,  Part  F). 

Concluding  Remarks 

A  method  for  the  analysis  of  transonic  flow  through  a  nearly  circular 
duct  with  only  gradual  charges  in  cross-section  has  been  derived.  One  is 
finally  lead  to  the  solutions  of  ordinary  linear  differential  equations. 
(Equations  (1.33)  and  (1.3^)).  These  solutions  can  be  carried  out 
numerically,  or  even  analytically  in  many  cases. 

Application  of  this  method  was  made  to  a  configuration  similar 
to  the  Wright-Patterson  Air  Force  Base  10-Jbot  Wind  Tunnel.  A  velocity 
distribution  through  the  duct  was  determined  and  the  disturbances 
caused  by  the  variation  in  duct  shapes  were  computed.  The  velocity 
profile  for  the  configuration  is  satisfactory.  The  subsonic  influence 
in  this  case  is  small  enough  that  the  usual  consideration  of  no  subsonic 
disturbances  propagated  into  the  supersonic  region  is  justified. 
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SECTION  IV 


APPENHX  A 


Derivation  of  the  Potential  Equation 

It  follows  from  the  vorticity  theorems  that,  with  certain  condition 
for  the  initial  flow  which  are  fulfilled  in  most  practical  cases,  there 
exists  a  velocity  potential  <p  (X,  7,  S),  where  the  coordinates  are 
shown  in  Figure  1.  (That  is,  the  velocity  vector  can  be  expressed  as 
the  gradient  of  a  scalar  function  $  ).  If  one  chooses  an  arbitrary 
function,  <p  (5t,  7 >  •)  and  determines  the  velocity  components  accordingly; 
i  .e  • , 

(3.1)  V*  =  <t>j 

(3.2)  V7=< f>? 

(3.3)  V  *  fes. 

"  7 

and  furthermore,  determines  the  pressure  from  Bernoulli's  equation,  then 
Euler's  equations  of  motion  are  automatically  fulfilled.  Therefore,  if 
a  velocity  potential  exists,  we  have  to  check  only  the  condition  of 
continuity.  The  continuity  equation  is  given  by 

If  one  inserts  Equations  (5.1),  (5.2),  and  (5.5)  into  Equation  (5.U) 


(}-5>  040-t- 


o 


The  derivatives  of  p  are  found  from  Bernoulli's  equation  for  three- 
dimensional  flow 


(3.6)  i.  +  -£-  (<pl  *■  ($>-*  + 


Con 


where  i  denotes  the  enthalpy,  so  that 

.  _ 

*■  ~  *+» 

From  the  isentropic  flow  relations,  since 

jL  -  /  £_  \X_' 

V?.  ) 

Consequently,  one  writes 


i.  =CyT for  perfect  gases. 


=  Os-O 

Or 

(5.7)  d  1  -  Q*  — -S- 
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If  Equation  (3.6)  is  differentiated,  and  Equation  (3.7)  is  inserted 
the  resulting  expression  is 

a1  + 

Or 

(3.g)  +■  4*^  ^ 4? +  4C^r)=° 

Writing  Equation  (3.8)  in  detail,  one  obtains 


(3.8b)  z  -  £-  (  d>_  (+)--  +  <f>-  0..  +  -V  (t^ 

>>y  ^yy  y*  y* 


(5.8c)  « -  -gi c»«  *».. »  4>?  <b ) 

Now,  if  the  partial  derivatives  of  ^  from  Equations  (3.8a),  (3.8b), 
and  (3.8c)  are  inserted  in  Equation  (3*5)>  one  obtains  the  expression 

0,9)  ~  I?  O*  <t>i«  +  <J>?  <l>Ry+4!!!-f^-)<l>*+?'t'« 

-  h  C*n  *  <to  <h?  +  %  +  1 


y  7 


+ -  f.  (^4^  * 


Lv  0?***  4>u>  WW^Q 


V.  y  y  7^  /  ^  ?*•" 

If  one  divides  Equation  (3*9)  by  ^  and  collects  the  terms,  one 
obtains  the  final  equation 

(3.10) 

OL* 

+ ^t-o-  -*  ftife*  -z  ^44^  ^-^,0 

This  is  the  complete  potential  equation  for  three-dimensional  flow. 
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APPENDIX  B 


Transonic  Law  of  Similarity 

The  general  three-dimensional  equation  for  a  steady  potential  flow 
is  given  by: 


(3 


,u)  (•-  +  f <.■*('-?%) 

-p  <t>u»  <b g u>  _  o  <pyto  $i: _ <fru>  <bv  _q 

y*  a>  Y  y’ci1  " 


The  velocity  of  sound 

(3.12) 


Cl  is  obtained  from  Bernoulli's  equation. 


<v>=  5±i 


y  | 

fc  -  T  W 


Where  the  total  stream  velocity* 


Let  us  consider  a  family  of  flow  patterns,  where  each  flow  pattern 
is  characterized  by  a  specific  value  of  a  parameter  z.  .  Let  the  velocity 
potential  be  given  by; 


(3.13) 

where 


4>  —  of  +  *  <9  ( y ,  **>)  +  z=x  Q  (x.y,  w)  + 


) 


(3.14)  *  =  % 

(3.15)  V  2:* 


This  hypothesis  is  justified,  by  the  fact  that  its  introduction  into 
the  differential  equation  for  $  ,  with  the  limiting  process  *-»0  , 

results  into  a  non-trivial  equation  for  •  Indeed  substitution  of 
Equations  (3.12),  (3«13)>  (3»14)>  and  (3*15)  into  the  potential  equation 
(Equation  (3.11))  gives 


(3.16)  <X* i:  (\-  ^ 


*-y-of  £  Cpt 


4 


* 

Ow 


Y 


c* 


ix 


The  value  of  the  sonic  velocity  as  given  in  Equation  (3.12)  dan 
be  written  in  terms  of  and  the  velocity  potential ,  and  is  finally 
expressed  as 


(3.17)  Ol*  —  d**  -  (  %-\)C^  <??)  o? 
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1. 


The  x  coordinates  of  corresponding  points  are  the  same 

2.  The  y  coordinates  of  corresponding  points  are  proportional 
to  2 

3.  The  difference  between  the  x  velocity  and  the  sonic  velocity 
is  proportional  to  £ 

The  difference  between  the  pressure  and  the  sonic  pressure  is 
proportional  to  * 

5.  The  deviation  of  the  streamlines  from  lines  of  constant  7  and 
<*>  is  proportional  to  i  t 

APPENDIX  C 


Boundary  Condition 

Along  a  bocfcr  at  rest  the  velocity  vector  is  parallel  to  the  surface? 
that  is,  the  component  of  flow  normal  to  the  surface  is  zero.  To  express 
the  relation,  let  the  unit  vectors  in  the  direction  of  increasing  x,  y 
and  w  coordinates  be  Ux,  u  u»  >  respectively.  The  total 

velocity  vector  then  is  expressed  by 

(3.2U)  <P%)  U.*  +  +  ^to*] 

The  vector  component  normal  to  the  duct  surface  is  given  by 

(3.25)  -  57  U-x  -  U-.  +U.y  =  0 

The  scalar  product  of  these  two  terms  is  the  component  of  the  velo¬ 
city  vector  normal  to  the  surface.  According  to  the  boundary  condition, 
this  component  is  zero.  Then  one  obtains 

-*  oT'  CPuj  —  Q 


(3.26)  -  0*0  *  ^  ^  +  O.*  “ 


The  assumption  that  is  a  small  quantity  implies  that  the 

slope  of  the  duct  and  surface  is  gradual,  and  consequently  ^  a 
are  smlUL  cjuantikios# 

Therefore,  the  expressions  and  <^<u,  are  of 

second  order,  and  can  be  neglected  in  Equation  (3.26).  This  results  in 


Ir. 


(3.27) 
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Derivation  of  the  Basic  How  Velocity 


The  basic  flow  velocity  is  the  velocity  obtained  from 

dimensional  flow  theory,  and  is  expressed  as  a  function  of  the  duct  radius 
at  the  cross  section  under  consideration. 

The  continuity  equation  is  written 

(3.2S)  r%  TT  p  V  =  r*  P  CX' 

where  V  is  the  absolute  value  of  the  velocity.  Since  ?V  is  a  function 
of  V,  it  can  be  developed  with  respect  to  (V  -  a*).  One  then  obtains 


(>.»>  =  --fc.oi.0 

If  one  introduces  the  relation 

(3.30)  Ar  -  v*-  v'q 

the  continuity  equation  is  rewritten  as 

]-n-(r.+  Avf  =  e  «*• 

If  ^  r  is  small^*  the  higher  order  terms  can  be  neglected  so  that 

(3.31)  Cv-<^  )  -  -qrj 

Since  ,  the  basic  flow  velocity  is  expressed  by 


(3-52>  CPb^Lw  n\T 


6--.  f 


APPENDIX  E 


Analytical  Representation  of  the  Constants  X. 


and 


An  analytical  representation  is  sought  for  the  integrals  that  occur 
in  the  expressions  X  0>*  an<*  YC.^*  Th®  ilrst  term  X.^is 

(3.33)  T.  =  C’v  £.^Cy)<ly  -  C‘ y  4y 

z 

Integrating^  yr  1 

(3.3M  y)  YUro 


5lbid.,  p.  11*5 
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Then^  since 


j.,  c-)„^  y )=t.  y) 

‘5-55)  iTsiU‘(%.v)+Vtvr)J 

The  values  for  the  arguments  of  the  Bessel  function  were 

chosen  so  that  ^ *  0.  However,  it  is  also  true  that  ^ 
so  that  Equation  (5.35)  simplifies  to  the  analytical  expression  ^**1 

<*•*>  Vt^%r) 

The  term  Koj*  is  defined  by 

\ 

(j.37)  k„k  =  £'  ys  2..^Cy)  <1/  =  y 3  X. (v  y)  4r 

If  one  denotes  the  argument  by  “•=<Vy>  Equation  (3.37)  is  rewrit¬ 
ten  as 


o  »< 


or  singly 


•58)  K.r=  'S.C-V* 


'°r 

An  integration^  gives  the  relationship^ 

«•»>  K„  =  w— u* 

or 

(3.m»  K o^=  5j7* J. (V) -  4 w  T,  (*)] ^'r 

If  the  Equation  (3.40)  is  rewritten  in  terms  of  y,  and  the  limits 
of  y  •  0  and  y  ■  1  are  inserted,  the  new  relation  becomes 

(5-w)  K.r=  L*  iC "3.  <1.,.  y) +  On.* - *  Y.r)7. (.%. r)] 

Since  —  =  O  the  K./.  5s  reduced  to 

0-w>  K.,.* 

^  Ibid.,  p.  11*5 
5  ibid.,  p.  145 
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The  term  X  is  defined  as 

(M3)  v  v)<iy 

^  / 

It  is  seen  that  the  integration0  leads  to  the  expression 

(3JA>  xvr-  ^  tv  y)-^v.,  (%.  y)X>,  (^,ry)]^ 

Also? 

(3.U5)  (■/■)  =  +  y-  i/y') 

(3.1.6)  -  2?+<  (y)  ^  2V Cy> 

The  substitution  of  Equations  (5.U5)  and  (5.2+6)  into  Equation  (3.2+2+) 
gives  the  new  relation 

(5.m)  T^=  £  tVC=)vry)+  "Sv  <^y\]. 

where  ~Sj  C^|v»a^  denotes  the  firwt  derivative  of  ^  v  with  respect 
to  the  argument. 

Since  the  arguments  are  chosen  such  that  3" v  (  vj*  ~  °  and 

the  integration  limits  are  y  *  0  and  y  ■  1,  the  value  of  Iv^  is 

(3.1^) 

The  term  K  vja.  is  written 

(3-1.9)  Kvr= jY+‘  £^CV>  Ay  =  £  yv+’  (y vr  y)  -Ay 

Substituting  the  following  differential  equation 

(3*5°)  J^Cy  Ty(;!ivry')]',-[y;:jvr-:Vj  ^yCXl^r'/)-° 

into  Equation  (3»l+9)>  one  obtains 

(3.51)  Kvr -  {■  y "  ^7  [y  V  y]} Jy + vt 

6  Ibid.,  p.  U+6 

7  Ibid.,  p.  12+5 


J.VJ3v(v)<«r 
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Integrating  by  parts,  one  can  rewrite  Equation  0.51) 

(j.52)  [-yvy  X ( V  ^  L  <*/] 

Further  integration  and  the  substitution  of  the  limits  of  y  ■  0 
and  y  ■  1,  leads  to  the  relation 

(5.53)  K,^=  ip  ("v  y)+  v  Tv  y)] 

Since  "3V  ('*')  the  anal3rtica-L  representation  of  V\v<^.  is 

(3.51.)  Kv>*  =  <£v  (%  V) 


APPENDIX  F 

Solution  of  Equation  (1.35)  if  ^  n«  is  Given  by  X  * 

The  solution  of  the  differential  equation  (1.55)  is  determined 
for  a  more  general  case.  Let  the  coefficients  of  the  differential  equation 
be  such  that  it  is  a  power  t  of  the  longitudinal  position. 

(5.55) 

This,  when  inserted  into  the  differential  equation  gives 


=.  O 


(3.56)  a*  %  i.  <  +  (*“•)  =  o 

In  order  to  remove  the  inconvenient  power  5^  in  the  coefficient 
of  a.*.  ,  let  us  introduce  a  new  independent  variable  u.  by  the 
equation  u.  =  *k  where  k  is  a  suitable  constant  which  is  chosen  later. 
Then  one  obtains 


(3.57) 


,  * 

Aye 

* 


-  k  XK' 


(5.58)  ^^.=Vvyt£k"° 

d*  _  x 

If  one  writes  Xs  and  inserts  Equations  (3-57)  and  (3 -5S) 

into  Equation  (5.56),  the  following  equation  will  result 


(5.59) 

du> 


ts.  +  r-S-  — 


x  * 
X  ^  %  1 


kV'^x) 
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where 


(5.60)  Y  = 

If  K  =.  \  -  ,  Equation  0.59)  simplifies  to 

•  a*'*-  *  A  *  X  * 

(5.61)  4  -v  X-  3  -v-  cl^-o 

du  ^  d*^  K 

Equation  (5.61)  has  the  following  form 

(5.62)  <X"  4-  (  +  (V  +  )  °^“  ° 

g 

the  solution  of  which  is  given  by 

(5.65)  G.=  US^(^) 

where  2  y  is  a  Bessel  function  of  p  order, 
let  us  set  the  following  relations: 

(5.6U)  ^  ~  2. 

(5.65)  p-S='-^p 

(5.66)  <?  =  (^i)1 

The  solution  of  the  differential  equation  (5.61)  involving  coeffi¬ 
cients  that  are  powers  of  x  is  expressed  as 

(5.67) 
where 

x  1-* 

(5.68)  S-— 


®  Ibid.,  p.  1^6 
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FIG- I—  CYLINDRICAL  FIG. 2-  SPACE  ELEMENT  IN 

COORDINATE  SYSTEM.  CYLINDRICAL  COORDINATE  SYSTEM. 


FIG.  3-  PROPAGATION  OF  DISTURBANCES  FROM  SUPERSONIC  JUNCTION- 
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FIG-4—  AXIALLY—  SYMMETRIC  DUCT  CONTOUR 


AXIALLY-  SYMMETRIC-DUCT-  CONTOUR 
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Figure  S.  -  Additional  Mach  Number  Along  the  X-Axis 


